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SUMMARY 


This  report  documents  the  method  which  was  developed  during  the 
creation  of  a  computer  program  to  calculate  thin  symmetrical  airfoil 
pressure  and  turbulent  friction  drag. 

The  computer  code  was  written  by  the  NAVAIRDEVCEN  to  provide  a 
consistent  method  for  the  determination  of  empennage  form  drag  of  high 
performance  aircraft. 

The  program  and  supporting  subroutines  are  written  in  FORTRAN  for 
use  on  a  CDC  6600  computer.  The  appendices  of  this  report  include  a 
sample  problem  input/output  listings,  computer  program  listing  and  a 
program  user's  guide  of  the  method. 
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INTRODUCTION 


This  report  describes  a  two-dimensional  linearized  flow  computer 
code  which  was  derived  by  the  NAVAIRDEVCEN  and  successfully  utilized  to 
calculate  form  drag  estimates  of  thin  symmetrical  non-lifting  airfoil 
sections.  The  computer  code  was  developed  during  a  NAVAIRSYSCOM  sponsored 
Aft  Fuselage  Twin  Nozzle  Fighter-Type  Aircraft  Drag  Project.  The  project 
involved  the  compilation  of  a  large  data  base  of  aircraft  geometric  and 
flight  performance  variations  associated  with  aft  fuselage  drag  measure¬ 
ments  of  various  wind  tunnel  fighter  models . 

One  of  the  problems  which  was  encountered  in  comparing  the  test  data 
from  the  various  sources  was  that  some  of  the  data  included  empennage  drag 
contributions  and  some  of  the  data  included  only  the  metric  tail  contribution 
of  interference  drag.  In  order  to  compare  the  data  on  a  consistent  basis, 
this  computer  method  was  developed  to  calculate  form  drag  estimates  of  the 
horizontal  and  vertical  tail  surfaces  for  the  wind  tunnel  models  with  metric 
empennage  sections. 

In  general,  the  analysis  can  be  extended  to  calculate  two-dimensional 
linearized  flow  pressure  drag  and  compressible  fully  developed  turbulent 
friction  drag  of  any  thin  symmetrical  non-lifting  airfoil  section  typical 
on  high  performance  fighter  aircraft. 


A  sample  case  is  included  in  Appendix  A,  together  with  the  resulting 
sample  output.  A  listing  of  the  main  FORTRAN  source  program  and  the 
required  subroutines  are  included  in  Appendix  B. 


DISCUSSION 


This  report  describes  a  two-dimensional  linearized  flow  computer 
method  for  pressure  drag  calculation  of  thin  symmetrical  non-lifting 
airfoil  sections.  The  method  approximates  the  drag  characteristics 
using  a  symmetrical  circular-arc  airfoil  in  lieu  of  the  actual  airfoil  section. 

The  transonic  pressure  drag  coefficients  of  two-dimensional  thin 
symmetric  airfoil  sections  can  be  represented  in  the  form  of  similarity 
parameters  (see  references  (a)  and  (b))  which  coalesce  geometry  variations 
into  a  unified  function  of  characteristic  similarity  coordinates.  The 
results  of  Spreiter  (reference  (a))  for  the  theoretical  and  experimental 
pressure  drag  of  complete  circular-arc  airfoil  sections  are  utilized  as  the 
characteristic  similarity  coordinate  within  the  transonic  linearized  flow 
regime . 
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The  transonic  drag  characteristics  are  extended  into  the  subsonic 
linearized  flow  range  by  matching  the  Prandlt-Glauert  drag  expression 
(see  reference  c))  to  the  low  range  transonic  similarity  function 
The  characteristic  transonic  similarity  function  has  also  been  extended 
into  the  supersonic  drag  range  by  matching  the  two-dimensional  symmetric 
circular-arc  simple  wave  drag  theory  (see  reference  (d))  to  the  high 
transonic  similarity  function. 

The  matching  in  the  above  analysis  yields  a  single  characteristic 
function  of  pressure  drag  for  thin  two-dimensional  symmetric  circular- 
arc  airfoil  sections  valid  throughout  the  tri-sonic  range.  The  character¬ 
istic  function  is  presented  in  Figure  1,  where  the  ordinate  represents 
the  similarity  drag  coefficient  and  the  abscissa  is  represented  by 
Spreiter's  similarity  coordinate.  The  utility  of  the  characteristic  function 
lies  in  the  ability  to  compress  the  pressure  drag  results  of  various 
circular-arc  profile  configurations  into  a  single  function  represented  by 
similarity  coordinates. 

The  tri-sonic  circular-arc  pressure  drag  characteristic  has  been 
programmed  in  terms  of  the  similarity  parameters.  The  computer  program 
presented  in  this  report  interpolates  the  characteristic  circular-arc 
pressure  drag  function  using  a  mathematical  spline  computer  routine  for 
various  input  data  values  to  yield  the  pressure  drag  coefficients. 

A  friction  drag  calculation  is  also  performed  within  the  computer 
program  for  fully  developed  turbulent  flow  across  a  flat  plate  flow  model 
which  approximates  the  thin  airfoil  profile  surface.  The  friction  drag 
coefficient  is  calculated  using  the  compressible  form  of  the  Prandtl- 
Schlichting  turbulent  flat  plate  equation  (see  reference  (e)). 

The  pressure  and  friction  drag  coefficients  are  combined  into  a  total 
form  drag  coefficient  which  can  be  referenced  from  the  section  airfoil 
area  to  any  other  reference  (i.e.,  wing  reference  area).  The  governing 
equations  that  were  utilized  in  formulating  the  pressure  and  friction  drag 
analysis  are  presented  in  the  methodology  section. 

The  assumptions  which  are  inherent  in  the  computer  program  associated 
with  the  derived  form  drag  methodology  can  be  summarized  as  follows: 

1.  Thin  symmetrical  two-dimensional  airfoil  pressure  drag  characteris¬ 
tics  can  be  approximated  with  sufficient  accuracy  by  an  equivalent 
thin  symmetrical  circular  arc  airfoil. 

2.  The  thin  smmetrical  circular-arc  airfoil  transonic  pressure  drag 
characteristics  can  be  represented  in  similarity  parameter  coordinates. 

3.  The  thin  symmetrical  circular-arc  airfoil  subsonic  pressure  drag 
coefficients  can  be  represented  by  the  Prandtl-Glauert  expression 

in  terras  of  similarity  parameters  matched  to  the  transonic  similarity 
coordinates. 
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FIGURE  1.  CIRCULAR-ARC  AIRFOIL  FRESSURE  DRAG  SIMILARITY  COORDINATES 
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4.  The  thin  symmetrical  circular-arc  airfoil  supersonic  pressure 
drag  coefficients  can  be  represented  by  a  wave  drag  expression 
in  terms  of  similarity  parameters  matched  to  the  transonic 
similarity  coordinates. 

5.  The  thin  symmetrical  airfoil  fully  developed  compressible  turbulent 
friction  drag  characteristics  can  be  approximated  with  sufficient 
accuracy  by  flat  plate  theory. 

The  computer  method  was  applied  during  an  Aft  Fuselage  Drag  Prediction 
Analysis  (reference  (f))  to  form  the  empennage  drag  estimates  for  several 
fighter  aircraft  configurations.  The  accuracy  of  the  method  for  the  thin 
symmetrical  airfoils,  from  the  above  analysis,  resulted  in  an  error  of  less 
than  ten  drag  counts.  The  demonstrated  degree  of  precision  is  considered 
acceptable  for  most  drag  estimates  of  high  performance  aircraft. 

METHODOLOGY 


Biconvex  Circular-Arc  Airfoils 


The  thin  airfoil  section  considered  in  the  analysis  of  this  report  can 
be  represented  to  sufficient  accuracy  by  the  following  equation, 

t32T[H-t)2]  <» 

The  theoretical  and  experimental  pressure  distribution  on  this  type 
of  airfoil  surface  has  been  investigated  by  Spreiter  (reference  (a))  for 
Mach  numbers  near  unity  and  for  various  thickness  ratios. 


Pressure  Drag  Coefficients 


The  pressure  drag  coefficients  are  obtained  from  the  known  pressure 
distribution  across  the  airfoil  profile  by  the  following  equation, 

5  2fi  tWt  <(-t)  <2 

•'0 

The  following  similarity  transformation  is  defined  for  the  transonic 
pressure  drag  coefficient. 


% 

c 


cd 


(3) 


The  transformed  drag  coefficient  has  been  experimentally  related  to 
airfoil  geometry  and  stream  flow  conditions  by  the  following  functional 
relation. 


"V  _ 

c  = 


f 


Moo  -  1 


t2/3  [(T+  1)  M.2]2/: 


(4) 


The  above  function  has  been  defined  by  Spreiter  for  thin  two-dimensional 
symmetric  circular-arc  airfoil  sections  valid  throughout  the  transonic  drag 
range  and  is  presented  in  figure  1. 
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The  subsonic  pressure  drag  coefficient  can  be  matched  with  the  transonic 
coefficient.  Taking  the  Prandtl-Glauert  pressure  drag  expression  (see  reference 
(b))  in  the  following  form, 

cd  »  0.0005  (5) 

Vl  “  M®2 

This  coefficient  is  subjected  to  the  same  similarity  transformation, 

^d  -  T(Y+  1)  _  0.0005 

=  7l - -  (6) 

t5/3  Vl  -  Moo 

The  transformed  subsonic  drag  coefficient  is  presented  in  the  labeled 
subsonic  portion  of  figure  1. 


Likewise  the  transformed  transonic  pressure  drag  coefficient  can  be 
matched  with  a  supersonic  drag  coefficient.  Using  the  following  expression 
for  supersonic  wave  drag  (references  (b)  and  (c))  for  thin  symmetric 
circular-arc  airfoil  sections. 


cd  = 


(7) 


The  transformed  supersonic  pressure  drag  coefficient  is  obtained  in 
the  following  form. 


cd  -  [(V  +  1)  Mg,2*]  ^3  # 

T5* 


(8) 


This  transformed  drag  coefficient  is  presented  in  the  labeled  supersonic 
portion  of  Figure  1. 
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Friction  Drag  Coefficients 

Fully  developed  turbulent  incompressible  friction  drag  of  a  thin 
circular-arc  airfoil  can  be  approximated  with  sufficient  accuracy  by 
the  flat  plate  flow  model  of  Prandtl-Schlichting  (see  reference (d) ) 
represented  by  the  following  equation, 

Cfl  ■  Uo^^-58  «0) 

The  compressibility  effects  can  be  accounted  for  by  use  of  the  following 
experimentally  developed  relationship  (see  reference  (e)), 


!i-  «  [1  +  .18m£  ]  ~*35 

C*i 


(11) 


Combining  the  compressibility  equation  with  equation  (10)  produces 
the  Prandtl-Schlichting  fully  developed  compressible  turbulent  friction 
drag  equation  in  the  following  form, 


cf 


[1  +  .18m2  ]  --35 


.455 


(log 


10  Re> 


2.58 


(12) 


The  thin  symmetrical  circular-arc  airfoil  pressure  drag  coefficient 
was  combined  with  the  compressible  turbulent  friction  drag  coefficient  to 
establish  a  total  form  drag  coefficient  based  on  the  airfoil  profile 
surface  area. 


CONCLUSIONS 

1 

The  computer  program  which  is  described  in  this  report  can  be  applied 
to  analyze  the  drag  characteristics  of  any  thin  symmetrical  non-lifting 
airfoil  section.  The  empennage  of  high  performance  aircraft  represents 
a  typical  area  for  analysis  using  the  method  of  this  report. 

The  computer  method  has  been  used  with  success,  for  several  years, 
at  the  NAVAIRDEVCEN  to  formulate  the  drag  characteristics  of  horizontal 
and  vertical  aircraft  tail  surfaces. 
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DESCRIPTION 

The  sample  case  presented  in  this  section  Is  representative  of  the 
engineering  problems  that  can  be  analyzed  using  the  thin  airfoil  form 
drag  method.  The  example  represents  the  vertical  and  horizontal 
airfoil  sections  from  an  F-17  empennage.  The  inputs  include  the  free- 
s t re am  mach  number  and  specific  heat  ratio  together  with  the  airfoil 
reynolds  number,  thickness  to  chord  ratio  and  the  ratio  of  airfoil  area 
to  wing  area. 

The  output  values  include  the  calculated  form  pressure  drag,  the 
turbulent  friction  drag  and  the  total  form  drag  referenced  to  the  wing  area. 
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APPENDIX  B 
COMPUTER  LISTING 


The  listing  provided  herein  includes  the  main  program,  subroutine 
spline  and  subroutine  LDAY,  which  are  currently  operational  on  the 
NAVAIRDEVCEN  CDC  6600  computer  system. 
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OPERATING  INSTRUCTIONS 


The  program  input  data  requirements  are  specified  in  Table  1,  where 
a  minimum  of  three  data  cards  are  required  for  execution.  The  first  data 
card  is  reserved  for  input  of  an  alphanumeric  case  identification  title 
in  columns  one  through  eighty. 

The  second  card  contains  the  airfoil  geometric  and  flight  parameters 
in  the  following  format;  the  free-stream  Mach  number  (XM)  is  input  in 
the  first  ten  columns,  the  specific  heat  ratio  (GAM)  is  input  in  columns 
11  through  20,  the  airfoil  Reynolds  number  (RE)  is  input  in  columns  21 
through  30,  the  airfoil  thickness  to  chord  ratio  (TC)  i»  input  in 
columns  31  through  40,  and  the  airfoil  surface  area  to  reference  area 
ratio  (ATAR)  is  input  in  columns  41  through  50.  The  data  supplied  on 
card  number  2  is  sufficient  for  the  program  to  calculate  the  required 
drag  coefficient  values  associated  with  the  input  data.  Additional  input 
data  with  the  same  case  identification  title  can  be  supplied,  as  necessary, 
by  repeating  the  card  2  input  format  for  each  set  of  additional  input  data. 

The  case  input  data  can  be  terminated  using  the  third  type  of  input 
card  listed  in  Table  1  with  parameter  ISWITCH  set  equal  to  "END"  in 
columns  1  through  3,  or  a  new  case  with  associated  identification  title 
can  be  initiated  by  setting  ISWITCH  equal  to  999.  Additional  input  case 
data  can  then  be  entered  following  the  input  specifications  of  cards 
1  through  3. 

The  output  format  consists  of  the  following  information;  the  case 
title  is  printed  as  a  header  in  the  form  specified  on  input  card  1,  the 
five  input  parameters  from  card  2  are  then  printed  followed  by  the  calcu¬ 
lated  values  of  pressure  drag  and  friction  drag  coefficients  based  on  airfoil 
surface  area  followed  by  the  total  form  drag  coefficient  based  on  the 
drag  reference  area. 
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TABLE  1 


PROGRAM 

INPUT 

DATA  FORMAT 

Card 

Quantity 

Mode 

Columns 

Description 

1 

Title 

Alphanumeric 

1-80 

Case  Identification  Title 

2 

XM 

Real 

1-10 

Free-stream  Mach  number 

GAM 

Real 

11-20 

Ratio  of  specific  heats 

RE 

Real 

21-30 

Reynolds  number  (in  millions) 

TC 

Real 

31-40 

Airfoil  thickness  to  chord  ratio 

ATAR 

Real 

41-50 

Airfoil  surface  area  to  reference 
area  ratio 

Repeat  card  number  2  format,  as  necessary,  for  additional  input  data. 

3  ISWITCH  Alphanumeric  1-3  For  end  of  data;  ISW1TCH  -  END 

For  new  case;  ISWITCH  *  999,  the 
enter  card  1  and  2  format  for 
additional  data. 
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